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Abstract 

A large class of 3-dimensional integrable lattice spin models is constructed. The 
starting point is an invertible canonical mapping operator !Zi t 2,3 m the space of a triple 
Weyl algebra. 7£i,2,3 is derived postulating a current branching principle together with 
a Baxter Z-invariance. The tetrahedron equation for 72-1,2,3 follows without further 
calculation. If the Weyl parameter is taken to be a root of unity, 72-1,2,3 decomposes 

( f) 

into a matrix conjugation operator Ri.2,3 and a c-number functional mapping 1Z\ 2 3. 
The operator Ri,2.3 satisfies a modified tetrahedron equation (MTE) in which the 
"rapidities" are solutions of a classical integrable Hirota-type equations. R.1,2,3 can be 
represented in terms of the Bazhanov-Baxter Fermat curve cyclic functions, or alterna- 
tively in terms of Gauss functions. The paper summarizes several recent publications 
on the subject. 

AMS subject classes: 82B23, 70H06; keywords: 3D integrable models 
PACS: 05.45-a, 05.50+q 

Introduction 

The subject of constructing and solving two-dimensional integrable lattice models is a 
well-established branch of mathematical physics. It has boosted important developments 
like conformal quantum field theory, quantum groups and the quantum inverse scattering 
method. But the analogous three-dimensional integrable lattice models have received little 
attention: the few known cases are complicated and direct physical applications seem to 
be remote. Many tools which are powerful for dealing with two-dimensional lattice models 
seem to have little relevance for the higher-dimensional cases. However, the Yang-Baxter 
equation can easily be generalized to tetrahedron- and D-simplex equations |21 B] which 
still guarantee the existence of commuting families of transfer matrices, and concepts from 
the quantum inverse scattering method like L-operators can be generalized to dimensions 
D > 2 Since many statistical systems require a truly three-dimensional treatment, 
exploring the possibilities of constructing and solving 3D-integrable models should not be 
useless. 

The first solution to the tetrahedron equation (TE) was found by Zamolodchikov in 
1981 1 , and the partition function of this N = 2-state model was calculated by Baxter jl] 

1 e-mail: gehlen@th.physik.uni-bonn.de; 2 pakuliak@thsunl.jinr.ru; 3 sergeev@thsunl.jinr.ru 
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using only symmetry properties, some factorization and commutativity. Later, Bazhanov 
and Baxter jS] showed how to generalize the Zamolodchikov model to N > 2, using 
the interaction-round-a-cube language (IRC-ZBB-model). Then, Sergeev, Mangazeev and 
Stroganov found the vertex- formulation of the ZBB-model These models share the 
feature the ZBB-Boltzmann weights are not positive definite, but at least in the thermody- 
namic limit the partition function becomes real. Baxter and Forrester [7| have investigated 
whether the integrable ZBB-model may describe phase transitions. They modify the ZBB- 
model introducing a symmetry breaking parameter (destroying the integrability). Calcu- 
lating the corresponding order parameter by variational and numerical methods, they find 
good evidence that in the integrable case the model is just at criticality (reminiscent of 
what happens in the 2D Potts model). 

There have been several attempts to construct physically more interesting 3D inte- 
grable models, avoiding the too restrictive Zamolodchikov TE, using "Modified TEs" 
(HI E3 E3 QUI- The notion of a Modified TE and the way to construct an integrable 
model from its solution has been first suggested by Mangazeev and Stroganov in [S]. 

A different starting point to find a large class of integrabel 3D-models has been used 
by Sergeev |17M18| . First an automorphism of a triple Weyl algebra is constructed, which 
satisfies a Kirchhoff-like auxiliary linear problem and a Baxter Z-invariance principle. This 
invertible canonical automorphism implies the TE, so this approach avoids the need of a 
tedious proof of the TE j^. If the Weyl parameter is taken to be the N-th primitive root 
of unity, a very interesting mechanism noticed before in |231 1201 1221 1131 1171 134j becomes 
operative: The mapping operator becomes a matrix conjugation (which is the quantum 
operator of an euclidian time evolution, or the Boltzmann weight) which depends on c- 
number parameters, which satisfy a set of Hirota-type classical integrable equations of 
motion. So we have a quantum integrable system which is parameterized by solutions of 
a classical integrable system. 

The different solutions of the classical integrable equations define various 3D-integrable 
quantum models. This framework contains the ZBB-model which is obtained choosing the 
trivial classical solution. The general solution to the Hirota-type equations for given global 
boundary conditions can formally be written down in terms of high-genus theta functions, 
using the Fay-identity |2T | 128 ) 12 ^ 130] . 

What is the use of these results? It is clear that the construction of a versatile integrable 
3D-statistical model with nice physical properties and its analytic solution is not yet 
round the corner, but the mathematical structures which such a model should incorporate 
become more clear. Whether this class of models is sufficiently wide to describe 3D phase 
transitions, will be seen. Another application which can be made right now is to use these 
3D models in an asymmetric geometry and to reduce them to new 2D integrable models 
|24j . Since Fermat curves are involved, these will be useful relatives of the chiral Potts 
model and of the relativistic Toda model [33] • The large variety of parameterizations which 
the MTE allows (in contrast to the Yang-Baxter equation) should help to overcome the 
serious parameterization problems which have prevented progress to finding more analytic 
results for the chiral Potts model. 

We shall not try to touch all interesting aspects of this program. We concentrate on 
explaining the construction of the basic triple Weyl automorphism and how it gives rise 
to the modified tetrahedron equation (MTE) and the classical Hirota-type equations of 
motion. We focus on the local properties and do not discuss boundary conditions and large 
systems. For the reader's convenience, we sketch some details of the derivations even if 
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these are rather straightforward. The interested reader may consult references |171 1361 133] 
which give more details on recent results on 3D integrable models satisfying MTEs. 




M 



Figure 1: Two cubes of the basic lattice intersected by the auxiliary plane (shaded) in two different 
positions: first passing through roi, tr>2, tt>3 and second through roi, xv' 2 , IU3. The second position is 
obtained from the first by moving the auxiliary plane parallel through the vertex A. The Weyl variables 
XV i, tv'i live on the links of the basic lattice. 



1 The canonial invertible mapping 72-123 

The models to be discussed in this paper are defined on a 3-dimensional oriented lattice, 
which can be imagined as a (possibly distorted) cubic oriented lattice. We use the vertex 
language: the dynamic variables tt)j live on the oriented links i of the lattice and the 
Boltzmann weights are associated to the vertices. We assume that the lattice is non- 
degenerate so that at all vertices exactly three lines meet. The Boltzmann weights can be 
considered as operators mapping the three dynamic variables on the incoming links to the 
three variables on the outgoing links. Fig^shows a perspective view of two adjacent cubes 
of the lattice (to keep the picture simple, we do not show the other vertices and lines). 
Only the dynamic variables on the links surrounding the vertex A are indicated. The links 
(or edges) of the lattice arise from the intersection of two planes each, the vertices from the 
intersection of three planes. We introduce auxiliary planes (shaded in Fig^) and locate 
the dynamic variables at the vertices of the lines where the lattice planes intersect the 
auxiliary planes. The mapping of the variables to, onto the to • can be regarded as moving 
the auxiliary plane through the vertex A. This way we may interpret the 3-dimensional 
statistical model as the euclidian time evolution of the 2-dimensional lattice which arises 
when we move the auxiliary plane. 

The next step is to specify the mapping, or, equivalently, the Boltzmann weights. 
We assume that the dynamic variables rOj are elements of the ultralocal Weyl algebra 
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Figure 2: Right hand side: magnified view of the vicinity of point A in Fig0 Left: magnified view of the 
auxiliary plane in the vicinity of r03 showing the four sectors formed by the lines rt>3 rt)2 and t03 roi . : 
We imagine a current \cj>3) flowing out of the intersection point into the four 2-dimensional sectors with a 
distribution steered by the Weyl variable tt>3 and the coupling K3. 



rOj € 2U® , defined as the tensor product of A copies of Weyl pairs (A is the number of 
vertices in the auxiliary plane): 

2rr® A = { U j = 1 ® 1 ® ... ®^u^,® ... , w» = 1 ® 1 ® ... ®^,® ... I % = 1, . . . , a } (1) 

i — th i — th 

place place 

with UjWj = q^WjUi, q G C. (2) 

In order to determine the mapping uniquely, we postulate a linear current branching 
principle (a kind of Kirchhoff law), together with a Baxter Z-invariance. Consider FigOU 
On the right hand side we show again the vicinity of the vertex A of Fig^ on top and on 
the bottom there are the two auxiliary planes. The three links of the basic lattice which 
connect the tt>i and the tt^, intersect in A. We now postulate currents \4>i) flowing out 
of the vertices into the surrounding four sectors of the auxiliary plane. The distribution of 
the current flowing out of rOj into the sectors is steered by the value of the variable 
XOi and one coupling constant k% , as shown in the left hand picture of Fig|2j The vertex 
i sends the current \4>i) into the sector left to the arrows, into the sector between the 
two outgoing arrows it sends g 1//2 Uj|0j) , to the right of the arrows Wj|</>j) , and between 
the incoming arrows KjUjWj|(/)j) . The total current in an internal sector is required to 
be zero and the total current reaching an external sector should be independent of the 
internal structure. In FigElwe show the lower and upper auxiliary planes of the previous 
two Figures. According to the rules just stated, the total currents flowing into the internal 
sectors \(p a ) and \(ph) must both vanish. \4> a ) receives contributions from \4>i), \4>2) 
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Figure 3: The triangles on the lower and upper auxiliary plane. 



and |03) - Since both structures differ just by a shift of one line through one vertex, there 
is a unique correspondence between the respective sectors. So these rules give: 



a) 


= wa|0i) + \<f> 2 ) + g 1/2 u 3 |</> 3 } 


= 


} b) 


= + w 2 |0 2 ) = w 3 |0 3 ) 




>c) 


= \<&) = |^> + |0 3 ) 




d) 


= \<P' 3 ) + q^WM = q^u.l 


fa) 


>e) 


= q 1/2 u' 3 \(j)' 3 ) = Kimwi \<j)\) + 


q 1/2 u 2 \<p 2 ) 


/) 


= Kiuiwil^} + K3U3W3I03) = 


- K 2 U 2 W 2 \4> 2 


'<?> 


= wil^) = W 2 |0 2 ) + K3U3W3 


\<Ps) 


h) 


= <? 1/2 u'i|0i) + k 2 u 2 w 2 |0 2 ) + 


w 3 |0 3 ) = 



(3) 

These eight equations determine the and uniquely in terms of the Uj and Wj and 
vice versa (we give details of the calculation in the Appendix). The resulting rational 
transformation can be written as: 



where 



w 1 

Ai 
A 2 
A, 



w 2 A 3 , 
A 2 _1 w 3 - 1 



u x 1 u 3 



W9 



u 2 



Wo 



A 2 - 1 ur 1 , 



A x 1 u 3 , u 3 = u 2 Ai 



(4) 



1 11 —l —l 

/ ' U x Wi + K\ WiU 2 



K 2 



U , Wo 



K 3 -1 -1 «1 K 3 _ x _ 

+ — u x w 2 — — u 2 w 2 

K 2 q 1 / 2 K 2 

J 1/2 U3W3 1 + K 3 W 2 1 U 3 . 



(5) 



L3 _ W1W3 

In (j3J the order of the factors on the right hand side is irrelevant since e.g. Ai contains 
only operators which commute with both u 2 and U3, in particular, it contains no w 2 , W3. 

The six eqs.© are written with only three operators Aj because the mapping conserves 
three centers: 



WiWn 



w 2 wi; 



UoU 



3 U 2 



u 2 u 3 ; 



Wq 



,/ —1 / 



U x = U1W3 



(6) 
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This rational mapping conserves the Weyl structure, i.e. it is canonical: from (J2J) it follows 
that 

u>; = r/" - w>:. (7) 
as may be verified by explicit use of (J3J) and ©• 

Definition 1 For any rational function of the Ui, ... , W3, the relations J2J) define 
the invertible and canonical mapping 72-1,2,3 

(7^1,2,3 $) (ui,wi,u 2 , ... ,w 3 ) d = ^(u^wi,^, . . . ,w 3 ). (8) 
7^-1,2,3 *s an automorphism of 2JJ® A and conserves the three centers 

The derivation of @, © from (j3J) given in the Appendix does not use that Uj, w-,- are 
Weyl operators. However, the canonical property of 72i,2,3 is obtained for ultralocal Weyl 
dynamic variables only. The inverse transformation has a similar form and we summarize 
it giving A^ 1 in terms of the primed variables: 

a-1 «1 / /-l 1/2 K 3 1 i-l . , i-l 

Al = %U3 -Q u l w l + K 3 U 2 W X 

K2 ^2 

A^ 1 = u 2 w 3 - g~ 1/2 k 2 u' 2 w 2 + uiw 2 , (9) 

a-1 K 3 i-l I 1/2 K l i-l 1 1 i-l 

A 3 = — w 1 w 3 — q 1 — u 3 w 3 + Ki w 2 u 3 

There is a way dual to © to formulate the current branching principle as "Linear 
problem". This is more convenient if one is dealing with a large lattice with definite 
boundary conditions. The resulting mapping is the same 1(1)) . 0, we give this alternative 
derivation in the Appendix. We refer to |36| I33j for a detailed exposition of the linear 
problem equations for large cubic lattices. 

2 Weyl parameter q a root of unity: decomposition of 1Z 
into a matrix conjugation and a functional mapping. 

Up to now, q was in general position. In all the following we fix q to be a primitive root 
of unity: 

q = to d = f e 2ni/N , JVeZ, N > 2. (10) 

Then and are centers of the Weyl algebra and we can represent the canonical pair 
(uj, Wj) by its action on a cyclic basis as unitary N x N matrices multiplied by complex 
parameters Ui,W{. The algebra 2U® A will then be represented in a A^ A -dimensional product 
space, see Because of the ultra- locality, the Weyl operators (uj, Wj) £ tttj will be 
represented trivially in all spaces j 7^ i, acting non-trivially only in the space with label i, 
where (now omitting the index i): 

u = ux ; w = wz ; 

I a) = \a mod N) ; ( a \ a') = 5 a}a >; x | a) = | a) uj a ; z | a) = | a + 1) . (11) 

In this representation the centers are represented by numbers: 

uf = uf , wf = wf . (12) 

We shall see that for q a root of unity, the rather complicated looking transformation 
can be represented in a rather simple way. 
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( f\ 

2.1 The functional mapping 7Z X23 

For q a root of unity, not only the N-th powers of the Weyl variables are numbers, but 
also the N-th powers of the operators Aj will be numbers. This is because in calculating 
the iV-th powers the cross terms drop out, e.g. (au + bw) N = (au) N + (bw) N due to 
ElV u j = 0, etc. We get 



Af = u x N u^ + u x N w? + «f< Uz N 



N N N N 

a N K l -N —N , K 3 —N -N , K l K 3 —N -N /,r,\ 

2 ^2 K 2 



Af = W * w -N + u^u^ + k£wZ n ^ 



So @ implies an analogous purely functional mapping of the Weyl centers: 

>N » N N >N a-N N >N a-N —N 

/N a-N -N fN a-N N /N \N N 

u± = A 2 w 3 , u 2 = A x u 3 , u 3 = A x u 2 ■ 



(14) 



Later we shall need not only the mapping of the centers (|14|) but also the mapping of the 
Ui, Wi itself onto the u^, w[. For this we have to take N-th roots and get a freedom to 
choose phases. If we demand the number products W1W2, u 2 us, uiw^ 1 to be centers of 
this mapping, then there are three free phases. 

(f) 

Definition 2 The functional counterpart of the mapping TZ x>2 s is the mapping TC X23 , 
acting on the space of functions of the parameters Uj, Wj (j = 1, 2, 3) 

(j^i%3°<Pj ( u i> w vU 2 ,w 2 ,u 3 ,w 3 ) d = (f)(u' x ,w' x ,u 2 ,w 2 ,u' 3 ,w 3 ) , (15) 

where the primed variables are functions of the unprimed ones, defined via 

u x = u x , w x = Wj , etc., (lb) 

such that w[w' 2 = wiw 2 , u 2 u 3 = u 2 u^, , u^u^ -1 = u x w 3 l . The three free phases of 

(f) 

the N-th roots are extra discrete parameters oflZ X23 . 
In obvious generalization: 

Remark 1 To any rational automorphism of the ultra-local Weyl algebra 2H® A J|J) there is 
a rational mapping in the space of the N-th powers of the parameters of the representation: 

N N 

2.2 Matrix structure of the mapping TZ\ } 2,3 

We now look at the matrix form of ^1,2,3 in the representation (|11|) . We separate the 
numerical factors from the constant diagonal resp. cyclic raising matrices Xj, Zj. Un- 
derstanding that the matrices Xj, act trivially in the spaces with j 7^ i, from © we 
immediately get the following complicated looking A^ 3 x A^ 3 -matrix equations: 

, M -l K X u' x _ x K 3 U X W 3 _ x _ x 1/2 KiK 3 u' 1 Wz _ x _ x 
(X X ) = X 2 H Xi Z 2 x 2 z 2 z 3 , 

' K 2 U 2 K 2 U X W 2 K 2 U 2 W 2 
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W 2 Wl 1/2^2^3 -1 , K 3 
Zi = — ZiZ 2 Z 3 -CJ 7 — Z2X3Z3 H r x 3 , 



.1 „.. „.l 



(x 9 ) = — Xi — o> 7 Xi Z1X0 H zix 9 x, , 

, M -l ^2-1 1/2^2^3 _l _l . K 3 W 2 U 3 _ 1 _ 1 

(Z 2 ) = Zo -LJ' Z, X 3 Zo H z, z 2 x 3 , 

w 3 101103 U>itt>2 

/ ""2^3 -1 1/2 U, 1'"2 _l . K 2 Wl 
x 3 = 7X1 X2X 3 lo ' r x 1 ZiX 2 H 7-Z1, 

/ M _l KlUiw'z _j _! K 3 W 3 _! 1/2 KlK 3 Ml^ 3 _j _i . 

(Zol = x ix 2 Zo H -z 2 — oj 1 x ix 2 z 2 . (17) 

K2U2W3 K2W2 K2U2W2 

However, eas. (|17j) can be rewritten in a compact way because there exists a (up to a scalar 
factor) unique iV 3 x N 3 conjugation matrix R.1,2,3 > such that for i = 1, 2, 3 

X- = R-1,2,3 Xj R-^2,3 ) z i = R-1,2,3 Zj ^2,3 > ( 18 ) 

and there is a closed formula for the matrix elements of R.1,2,3 which later will be given 
explicitly. Let us write the matrix elements of Ri 2 3 in our representation (|11|) as 

<fl,i2,i3|R 1 ,2,3U 1 ,i2,i3) = 4lfJ:- (19) 

Then e.g. the from the equation for x 3 in (|17jl we see that we have to find Rll'j^'i^ such 
that it satisfies the recursion 

( , ,33 _ !f2, ,-11+12+13 1 T>h,h,33 _i_ ,-l/2+i 2 -ii^l^2 R ii,i2-l,i3 _ K 2^1 ^31,32,33-1 _ n 

V ui ^1,^3 +tJ ui^^'W u' 3 

(20) 

together with five other similar relations from (|17[). The indices are cyclic mod iV. 

Now, going back to we remark that in our representation also any rational function 
$ of the ui , . . . , w 3 is represented by a iV 3 x iV 3 matrix with entries which are rational 
functions of u%, ... ,11)3 . The mapping 7^2,3 transforms the matrix entries by 7Z\ v ; 3 
and conjugates the matrix by Ri,2,3 • So we get the 

Remark 2 At q a root of unity, the mapping 7Z\ t2 ,3 can be represented as the superposi- 
tion of the pure functional mapping TZ\ ^ 3 an d the finite dimensional similarity transfor- 
mation Rl,2,3 •' 

^1,2,3 o d> = R 1)2i3 (n[ f l 3 o $) R-i )3 . (21) 

This remarkable property arising in models where the quantum variables are elements of 
a Weyl algebra at root of unity, has been pointed out in |23j and in ^M,, |20| . 



3 Tetrahedron equation and modified tetrahedron equation 

We postpone for a moment the discussion of the explicit form of Ri,2,3 and of the possible 

( f) 

form of 1Z\ 23. We shall first see that, once we have constructed the invertible canonical 
mapping of the triple Weyl algebra 7^1,2,3, no calculation is needed to conclude that T^i,j,k 
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satisfies the Tetrahedron equation. We want to keep the point of view taken in Fig^that 
the mapping 7^1,2,3 comes about by moving the auxiliary plane through a vertex of the 
basic lattice. For considering the single mapping 7£i,2,3 we selected three planes of the 
basic lattice and intersected these by the auxiliary plane. Now we look at four (with no 
pair being parallel) lattice planes intersecting the auxiliary plane. In the auxiliary plane 
this gives rise to a figure called quadrangle shown in Fig^J A quadrangle consists of four 
lines which form two triangles and one four-sided figure. There are six crossing points. 
Shifting the auxiliary plane within the basic lattice, we can reverse (from A to V) either of 




A 



Figure 4: Quadrangle in the auxiliary plane formed by the directed intersection lines of 4 oriented lattice 
planes. 



the two triangles: Shifting such that the plane producing the line B gets to the other side 
of point 1 (this happens moving the auxiliary plane through the vertex formed by planes 
B, C, D) we perform the mapping 7£i,2,3- Alternatively, we may shift the plane C through 
the point 6 (moving the auxiliary plane through the vertex of A, B, C), so performing the 
mapping K 3:5fi . 

As shown in Fig 03 there are two different sequences of mappings which, starting from Qi, 
lead to the same quadrangle Q$: either Q\ —* Q2 — ► Q3 — > Qa — > Q5 or Qi — > Q$ — > Q-j — > 
Q& — ► Qb- Since the mapping T^i jk ls invertible and canonical, we conclude that the two 
corresponding products of the 1Z must be the same mapping: 

7^123 • 7^145 • 7^246 ' 7^356 ~ 7^356 • 7^246 ' 7^145 ■ 7^123 , (22) 

i.e. that Ttij k satisfies the Tetrahedron equation 2 . (|22|) is an operator equation in the 
space of the six Weyl variables. An equivalent interpretation is that the six dynamic 
variables live on the six edges of a tetrahedron, and the right hand side of l|22[) is obtained 
from the left hand side by moving one corner of the tetrahedron through the plane defined 
by the other three corners. Inserting the decomposition 1)21(1 of 7Zij,k into a matrix- and 

2 We use the following notation for the superposition of two mappings A and B: 

{(A-B)o^) = f (_4o(6o$)) . (23) 
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Figure 5: Graphical image of the two equivalent ways of transforming the four- line-graph ("quadrilat- 
eral") Qi into graph Q5, which leads to tetrahedron equation. Observe that each graph contains only 
two triangles which can be transformed by a mapping 1Z: In graph Qi either the line 124 can be moved 
downward through the point 3 (leading to graph Q2), or the line 456 can be moved upward through point 
3 (leading to Qs). Both the left hand and right hand sequences of four transformations lead to the same 
graph Qs- 



a functional part into (|22l) . we get 



R 



123 



^123 "! R 145 ( 7^ 



(/) 
115 



R-246 I ^ 



R 



n 



if) 

356 



R246 ( 



if) 
246 



1^356 °^ 



R 35e) } 



R 



-1 

246 



R 



145/ 



R 



-1 

123 



R-145 {^-145 ( R 123 (^-123 ° *) ^h) } R 145 ) R 24e} 



R 



-1 



356- 

(24) 



In (|24|) we have to apply the functional transformation repeatedly to matrices Rjj,fc and 
1 f) 

to KY/u and $. Let us write shorthand 



R« 

R( 4 ) 
R( g ) 



R 



123; 



R( 2 ) 



7?(/) 7?(/) 
/-t -l,2,3 /-t -l,4,5 /c 2,4,6 



~~ 'M,2,3 
R-356! 



R 



145; 



R( 3 ) = 
R( 5 ) = 

^5,6 R 246 5 



^•lfi^^M^ ° R 246; 

7?(/) -p(/) -r,(/) oR . 
'^3,5,6 /,d 2,4,6 'M.4,5 ° K 123. 



R< 8 > 



R 



356- 



(25) 
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Then l|24|) can be rewritten as 

= (r( 8 )r( 7 )r( 6 )r( 5 )) {ngln^-n{%-n{g^} (rWr^r^r^) -1 . (26) 



Since, as has been discussed in Sec. 12.11 any rational automorphism of the ultra-local 
Weyl algebra implies a rational mapping in the space of the N-th powers of the parame- 
ters of the representation, its is a direct consequence of dm) that the K^F) of (JTSJ) solve 
the tetrahedron equation with the variables Uj, Wj , j = 1, . . . , 6: 

-123 /c 145 '^"246 /c 356 /c 356 '^246 /c 145 '^123^ ^V^l ) • • • ) ^6 ) ~ u - 

Using Maple, by a ten-line-half-minute-program one can easily check the validity of this 
functional equation straight from the definitions (|14|) . I|13[l. 

Now the R J depend not on the u$ , «^ but on the Uj, Wi, see (|17|). So we have 
to chose the phases, which are at our disposal in defining TZ\ ■ k such that it solves the 
functional tetrahedron equation for rational functions of the variables Uj and Wj directly. 
However, not all phases of the Uj, Wj can be chosen independently because there are 
various centers which should be conserved. We shall not enter into the discussion of these 
phase choices here, in |35| it is shown that 16 phases can be chosen arbitrarily. Once a 
consistent choice of phases has been made, we obtain (|27|) in the variables Uj, Wj and can 
cancel the central functional terms in (|26j) . So we arrive at the 

Proposition 1 The finite dimensional matrices R satisfy the Modified Tetrahedron Equa- 
tion (MTE): 



R-1,2,3 ' ( ^-1^,3 ° ^1,4,5) • (^l{2,3^L4,5 ^2,4,6) ' \^,2,3^1,5^%l,G ° R 3,5,6 



R-3,5,6 • ( R 2,4,6j ' ('^■/,5,6 R %Ig R l,4,5j ' {^g}, 6^2^^1,1,5 R l,2,3 



(28) 



The left and right hand sides of (|28jl may differ by a scalar factor, which arises when we 
pass from the equivalence of the mappings to the equality of the matrices. So, writing all 
matrix indices and using again the abbreviations ([25)1 . the MTE reads 



£ 

31— h 



R{1) yuh,j 3 / (2) xfcij4,j5 r R{3) \ k2 ' k4 > j6 r R{A) \^MM 



11,12,13 v ' 3 1,14, «5 v ' 32,34,16 v ' 33,35,36 

33,35,36 / ,^\h,H,k(i / t( ,,\jl,k4,k 5 / , K ,\ kx ,k 2 ,k 3 

(29) 



31-36 



«3,«5,*6 x / «2,MJ6 v / UJ4,J5 v ' 31,32,33 



Here /? is the scalar density factor. Taking determinants in (|29|) we can express the iV 3 -t/i 
power of the scalar factor as 

^3 _ detR (1) detR (2) detR (3) detR (4) 
P ~ detR< 8 > detR< 7 ) detR( 6 ) detR( 5 )' { ' 

This can be obtained from the determinant of one single matrix Ri 2 3 just by substi- 
tuting the respective arguments. The components of the eight R^^ may be considered as 
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Boltzmann weights of the models. However, in general these matrix elements will not be 
positive. We call eqs. (j28j) or (|29|) modified tetrahedron equations because in (|28j) the Rw 
depend on several " rapidity" variables u\ , . . . , wq which are not the same on the left- and 

right hand sides of the equation, but rather are related by various functional mappings. 

(f) 

We shall later see that there is much freedom in choosing 7£^ 2 3> accordingly leading to 
different models. 



4 Details on the matrix- and the functional transformations 

4.1 Matrix part of 7^i 2 3 at root of unity in terms of Fermat curve 
cyclic functions W p (n) 

In Sec l2.2l we introduced the conjugation operator R-1,2,3 an d its representation -R^^fg > 
which allowed to write eas. (|17|) in the form ()18j) . We just have seen that the matrix 
elements of the Ri,2.3 play the role of Boltzmann weights and solve the MTE. 

In (I2UI) we had remarked that Rj}'^ 2 '] 3 has to satisfy six recursion relations with respect 
to its modulo N defined indices, but have not yet seen its explicit form. 

We now show that the Ftf 1 '? 2 '- 3 can be constructed from Fermat curve cyclic functions 
W p (n) which were introduced by Bazhanov and Baxter [3] for building up the Boltzmann 
weights of their model. The form (|33j) which we shall use has first been given in several 
years ago. 

Introduce a two component vector p = (x, y) which is restricted to the Fermat curve 
x N + y N = 1. Then define the function W p (n) by 



Wp(0) = 1, W p (n) = Yl - J for n > 0. (31) 

u=l 



W p (n) satisfies the obvious recursion relation 

W p {n - 1) = (1 - to n x) y- l W p {n). (32) 

and is cyclic in n : W p (n + N) = W p (n), because of n^o 1 0- ~ ^ x ) = U N ■ Observe 
that we need the Fermat curve restriction in order to get the cyclic property. In order to 
satisfy the six recursion relation required by (|TT|). -R^'^'/g 3 must be taken to be a kind of 
cross ratio of four functions W p , depending on four Fermat points p%, . . . ,p^: 

Proposition 2 In the basis ijiij) the mapping relations H8\) are solved by the matrix 
where the x- coordinates of the four Fermat curve points are connected by 

X\ X2 = UJ X3 X4 , (34) 

and the Fermat points are given in terms of the variables uj, Wj, Kj, j = 1,2,3 by 

-1/2 u 2 -l u 2 -1 K 2 U 2 

U 1 Ul Ki U 1 





uj- 1 / 2 


U-2 




XI 






, x 2 = 




Kl 


Ul 




m 


IL>1 




VA 




= Kl — 




— = UJ 


Vi 


u 3 




yi 



1/2 „ ^3 V3 W 2 y 4 1/2 ^3 U 3 

K3 — , — — — , — — to T , [6b) 

W2 2/2 ^3 m ki w 1 
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The u'j, w'j and Uj, Wj are related by the functional transformation \15\) : 



■3 - -j. w j = ^i2,3 ° ( 36 ) 

In order to prove this proposition, we have to check the six recursion relations implied by 
(|17J) and have to compare the Fermat point coefficients arising from (j32j) with the Weyl 
coefficients in (|17|). We give one example: 
satisfies 



t>3i,32,33 
l il ,12+1,13-1 



r>3l,32,33 
il,«2,«3 



2/1 

2/4 



1 



o/2-jl+l^ 



1 - ^2-il+l 



3'1 



which can be rewritten in the form 

Hi ,12 ,13 



•a+i, 



t>31,32,33 
il,l2,«3 



£i2/4 



nJU2,j3 
U,i 2 +1,«3-1 



2/4 



T>31,32,33 

«1,«2 + 1,«3-1 



X 4 0^1X4 "»X,»2-rx,»3-i W «a+ia; 4l/l n^t-x^-x ( 37 ) 

Using the identifications (|35|) we see that this agrees with the matrix elements of the first 
equation of (|17j). 

In order to find the density factor of the modified tetrahedron equation, in (|30|) we 
have seen that we need the determinant of Ri 2,3- Its calculation is made nontrivial by 
the presence of the factor 5i 2+ i 3 j 2+ j a u;^ 1 " 11 ^ 3 in (f33|) . There is an interesting technology 
involving the Fermat curve functions W p (n), details of which can be found in [51 111! 131)] . 
compare also the appearance of W p (n) in connection with the quantum dilogarithm e.g. 
in eq.(3.7) of [201 ■ The calculation of det R is given in detail in [3SJ ■ The result is: 



det R 



N 



X4 

.2/1 2/2 



V(x 3 )V(x 4 ) ) 



(38) 



with V(x) 



(i 



v+1 



a:)". 



4.2 The matrix i?f f f for N — 2. 

n,i2,«3 

For N = 2, using combined indices £ = 1 + *i + 2i2 + 4^; j 
give Rll'i^'ig explicitly in matrix form: We define 



1 + ji + 2j2 + 4j 3 we can 



Y 



and get 



Rf 









1-Xj 




for 


i = i, 


2, 3, 4; 






1 + X j \ 


1 1+Xj 








Yi 
Yj 


for ij 


= 13, 14 


23, 


24; 


Z12 = 


YxY 2 ; 


Z34 


/ 


1 


Z24 














Z23 


—Z34 




Z\3 


Z13Z24 














— Z\2 


Z14 










Z13Z24 


Zx3 




— Z\2 
















Z24 


1 


— Z34 


Z23 
















Z23 


Z34 


1 


— Z24 
















Z\2 


Z14 


— Z\3 


Z13Z24 










Z\A 


Z\2 














Z13Z24 


— Z13 


\ 


Z34 


Z23 














—Z2A 


1 



(39) 



(40) 
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The determinant can be calculated directly: 



det R = (Y\ Y 2 (1 + Kf 2 )(l + y 4 ~ 2 )) 4 . (41) 

4.3 The functional mapping in terms of trilinear Hirota equations 

When we wrote the MTE in SecEl we did not specify the details of the functional trans- 
formations TZ\ ■ k which act on the scalar parameters of the quantum operators Rij,fc (by 
(|33|) the scalar parameters then determine the Fermat variables). We are interested to 
write this more explicitly and we like to understand which restrictions the MTE imposes 
on the Tl\2k- 

Slightly rewritten, eqs. (|13j) . (|14|) which define the functional mapping, are 

/ N N N N N N 

U 3 U 2 _ , \ N _ U 2 W{ U 2 N W{ 

~ W ~ ~ TW ~ \ U 3 A l u 2) — ~ ^ N N + K l — W ' 
U3 U 2 u l u l u 3 u 3 

N N NNNNNNNN 

u i w 3 a ^^ r Kfuf ftj W3 «i «3 ^1 W3 



1 X , X ( w 3-h-2 u l) at at + \' y 1 \ \ \' 

u[ W 3 «2 U 2 K 2 w 2 K 2 u 2 w 2 

I N N N N N N 

W l W 2 / 1 A \N Wn WnUn AT Un 

If we just use (|l2]l repeatedly to calculate the in (|25|) . we get very lengthy expressions. 
R( 3 ) will be much more complicated than R^ 2 ) etc. It is not transparent, which variables 
are independent because the conservation of the centers is hidden in complicated formulas. 
However, the nice symmetry of Fig|S] suggests that a different parameterization should be 
possible which treats all operators in l)28|) on essentially equal footing. Indeed this is 
achieved applying a Legendre transformation or, equivalently, by the introduction of tau- 
functions as it is common practice in the theory of classical integrable systems Let 
us see in an elementary way how this works. 

In FigEl we show a magnified view of the two quadrangles Q\ and Q 2 of Fig03 
Consider Q\ : The twelve scalar variables Uj, Wi (i = 1,...,6) are located at the six 
vertices labeled i = 1,...,6. The vertices cut each directed line A, B, C, D into four 
line-sections. We label these line-sections as shown in the figure. Ratios of line sections 
will be our new variables and we explain the rule for the variables u\, w\ associated with 
vertex "1" in the upper left corner: 

• u\ = c 2 /c3, i.e. a Ui is expressed as the ratio of the sections before to after point "1" 
on the right pointing line through the vertex " 1" . 

• w\ = d^/d 2 , i.e. a Wi is expressed as the ratio of the sections after to before point 
" 1" on the left pointing line. 



This way the twelve scalar variables u\,wi, . . . ,uq,wq of each quadrangle are expressed 
in terms of eight "internal" line sections, 01, a 2 , . . . , d\, d 2 for Qi , and eight "external" 
line sections do, 03, . . . ,do, ofe . 

Passing from one quadrilateral to an adjacent one, corresponding to a mapping 1Z, 
always three of the "internal" lines are changed, and for distinction, we attach to these 
changed variables dashes (if left-turning in Fig|S] and daggers (if right turning). The 
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Ql B C B Q 2 




Figure 6: Parameterization of the arguments uf , w\ of the rational mapping TZ?f, ' & in terms of 
line-section ratios. 



"external" line sections remain unchanged in all mappings and so are no relevant dynamic 
variables. 

( f) 

Applied to the transformation 1Z\ 2 3 this means that we substitute 

c 2 b 2 h d 3 d 2 c 2 

u\ = — ; u 2 = t-; u 3 = — ; wi = — ; w 2 = 3-; w 3 = — ; 
c 3 63 b 2 d 2 di c\ 

1 c i r b i r b' 2 / d 2 , d 3 , C3 

Ul = ^ ; U2 = ^ ; U3 = ^ ; Wl = ^ ; W2 = 4 ; W3 = 7 2 > (43) 

and (|42[1 become tri-linear equations for the three dashed variables: 

b 2 > N cZd% = + + 

K 2 ^2 c 2 ^2 = K ib 3 d 2 + K 3 b 2 c% d^ + Ki K^b^ c 2 d^ ; 
b%c»d 2 > N = b»c»d» + b?c?dZ + ^bM4; (44) 

The merit of this parameterization is that it automatically incorporates the invariance 
of the centers of the mapping: w^w^ = wiw 2 , u 2 u' 3 = u 2 u 3 , u^w'^ 1 = u\w^ 1 . It is 
also simple to express the Fermat variables appearing in (|35j) in terms of the line-section 
parameters. For R.1,2,3 we label the Fermat coordinates by a superscript M (as in (|25|)). 
and get: 

(1) -1/2 b 2C3 (1) 

x] ' = oj ' — ; x\ ' = to 

Klb 3 C 2 

(!) L 7 (1) 

V 3 _ K l03Q3 j/4 
y CD b' 2 d 2 ' y CD 

Now to find all the functional mappings which appear in the MTE (|28|) . we have to write 
also the other seven relations corresponding to Q 2 — > Q3 , . • • , Qi — * Qsi which look all 



1/2 K 2 bic' 2 _ 


r (l) 


- o;~ 1&lC3 


6' 2 ci 




6' 2 c 2 


i/ 2 n 3 c 2 di _ 




cid 3 


Cl^2 




c 2 d' 2 ' 
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similar to (|44|). (|45|). E.g. Ri^s leading from Q 2 to Q3 is governed by a 2 , c'{, d'(, c' 2 , d 2 etc. 
which are related by 

n "NNrN _ NjNjN , „NjNj IN , „N n N „N j IN. 
a 2 c l % — a l c 2 "1 ~r a 2 c 2 ^2 + K l a 3 c j U 2 , 

K4 (2 2 = Ki (Z3 Cq 

K 5 a 2 c 2 ^5 ^3 "0 ' 

a%c?d" N = a%c$d 2 ' N + a?c$d? + ^a?c?do- (46) 
( f) 

Since we know already that A is a solution to the functional TE, one may easily show 
that eight of the 8x3 trilinear eqs. of type (|44|) are superfluous. For the whole MTE we 
have introduced 16 line sections to replace the 12 scalar variables, but the eight "external" 
variables ao, • • • ,^3 and the six couplings turn out to be irrelevant, i.e. these can be 
eliminated by a re-scaling |35j . 

So the conclusion is that the eight "internal" variables a\, a 2 , . . . , d\, d 2 together with 
sixteen phase choices parameterize the MTE eq. (|28j) or (|29|) . This is much more compli- 
cated than the rapidities which appear equally on both sides of a Yang-Baxter-equation. 
Of course, there is no difference property of the parameters of the MTE. 

4.4 Solving the MTE using tools of algebraic geometry 

It is well-known [23 1231 EH1 1301 131] that discrete integrable classical equations can be 
solved by methods of algebraic geometry. The bilinear Hirota equation is equivalent to 
Fay's identity for theta-functions on the Jacobian of an algebraic curve [30] • The tri-linear 
equations (|44f) are also of the correct type to be solved using Fay's identity. 

We use standard concepts, see e.g. [321: Let T g be a generic algebraic curve of genus 
g and to the canonical g-dimensional vector of the holomorphic differentials. For any two 
points X, Y G T g let I : T® 2 C 9 be the Jacobi transform: 



l Y x = I(X,Y) = f f to G Jac(r g ) 
Jx 



(47) 



Let further E(X, Y) = -E(Y, X) be the prime form on T® 2 , v G C 9 , and 6(v) the zero 

characteristic theta-function on Jac(r 9 ). Then, for A, B, C , D G T g Fay's identity 
states 

e(v) 6(v + ii + i£) = e(v + i£) e(v + 1^ 



C . , tA^,„. , rC,E{A,B)E(D,C) 



E(A,C) E(D, B) 

In order to solve (|44l) . (|46|) . etc. we can build a tri-linear identity by applying (|48|) twice: 
G(v 4 1%) 6(v + I? + if) 6(v + If + 

- U(v 4 l x 4 l z ) U(v 4 l y j U(v + l z + L Y ) 

4 G(v 4 I? 4 IV) 6(v + Ig) G(v 4 4 4 If) 

- 0(v 4 I x 4 I y + I z ) 0(v 4 I y ) 0(v + l z ) E{z , fX)E( Y, iX)E( Y tZ) • ( 49 ) 
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Q, X, Y, Y' , Z, Z' are arbitrary distinct points on T g . The dependence on Q is trivial since 
it appears only in a fixed combination with v. It is not easy to handle the prime forms 
themselves, but here only cross ratios of prime forms are needed and these are well-defined 
quasi-periodical functions on r® 4 : 

E(A,B) E{D,C) _ e eodd (i*)e eodd (i c D ) 

E(A,C)E(D,B) 6 £odd (lS)G eodd (lg)' 

where &e odd is a non-singular odd characteristic theta function such that £odd (O) = 0. 

Consider eas.(|44|) for showing the application of (|49|) to solve our tri-linear equations. 
If we identify 

6f ~ G(v + I Q X ); ~ 9(v + 1% + Iy"); 6f ~ 9(v + I Q X + IV + if); 
b' 2 N ~ G(v + I x + if); £ ~ 6(v + if); cf ~ 8(v + Ij + If); 
^ ~ G(v + if + Iy"'); ^ - 6(v + if); G(v + if + if). (50) 

then, not discussing here the prime form cross ratios, the first eq. of (|44|) takes just the 
form (|49j) . The third eq. of ([44 [I takes this form if we exchange X <-> Z in (|49|1. In order 
to obtain the second eq. of replace in (|4*9"|) X — > Y' , Y — > X and v — > v + Iy . 

Here we shall not give the details |38| of the calculation which shows that all 24 tri- 
linear equations can be written in the form (|49j). For parameterizing the line section 
parameters and so parameterizing the eight R^) which appear in the MTE, we have to 
choose a T g and to introduce just eight points X, X' , Y, Y' , Z, Z\ U; U' G T g . 

From these formulas simple parameterizations can be obtained taking the rational 
limit. Then the theta functions become trigonometic functions and the cross ratios of the 
prime forms become simple rational cross ratios |36j . 

5 Free bosonic realization of Ri 2,3 

Together with introducing the ZBB-model in |2H], Bazhanov and Baxter also considered 
a related continuum bosonic model. Since our framework includes the ZBB-model, we 
expect that such a bosonic representation should also exist for R-1,2,3- Since this broadens 
the class of models, we shall present here the basic formulas. We shall not consider this 
realization a limiting case of the discrete realizations, but rather a new representation. 

Instead of the cyclic weights W p (n) of 1)31)1 we introduce the following Gaussian weights: 

WJo) = exp f \ — — a 2 ) ; a G R, x G C; 3m — >0. (51) 

\2n x — 1 / x — 1 

At each vertex j of a graph define a pair of operators q^, pj satisfying [q^, pj ] = ififijji, 
and scalar variables Uj, Wj. We choose a basis \aj) with (o~f\(Tj) = <5(ay — crj), such 
that for tp(a) G I? we have 

M?,) = f fa#>; (o-il qj IV) = o-iV'(^); fol P. IV) = y (52) 

For each set of vertices we define the corresponding operators and scalars and a direct 
product space of the single vertex spaces. Consider now the following mapping R.1,2,3 in 
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the product space of three points j = 1,2,3 by 



_j dlogu'j dlogu'j _ / 

R-123 R-123 = / "31 ^fe + "51 Pfc = °l? 

^^l ogUfc <91ogw fc J 



fc=l 

i dlogw', <91ogw/- 

R-123 Pj R-123 = 2_«1 + ol Pfc = Pr ( 53 ) 

fe=l 

where the relation between the primed and unprimed scalar s Uj, wj is given by ()15[) with 
iV = 1 . The Kj are further parameters ("coupling constants") at the vertices j. As the 
following Proposition shows, equations (|53j) are the bosonic continuum analogs to eas. (|17|> 
of the discrete case. 

Proposition 3 In the basis \5 6 J\) the operator Ri 2,3 has the following kernel 
(o"i , °"2 > °3 fRi,2,3 ki , 02 > a 3) 

= 5{a 2 + a 3 -a 2 -a 3 ) e 1 1; s _____ -7-, 

H/ X3 (cj 2 - crj kV- 4 (cr 2 - cr-J 

u;ii/t i/ie constraint x\ x 2 = 14. In terms of the variables Uj, Wj, Kj, (j = 1, 2, 3), i/ie 
Xfc S K are defined by (obtained by putting formally cj -1 / 2 — » —1 in ^5j) ): 

1 U 2 U 2 U 2 K 2 «2 ,ce\ 
x l — i 2^2 — ~ K 2 T j ^3 — ) X4 — r j \oo) 

where u\ an d u 2 are defined as in \3b)) with N = 1. 

Proof : We give the proof for one of the six equations ()53|) , as the other equations follow 
analogously. Let us write shorthand |c) for \o\, a 2 , 0-3) and d 3 a = da±da 2 da 3 etc. We 
consider: 

Jd 3 a' (a\H 1>2t3 \a l )(a'\q. 3 \a") = I d 3 a' ( a | q£ | a' )( a' | R 1)2 , 3 | a") (56) 

which should be satisfied for all ax, o~2, o~ 3 , o~'{, o 2 , a 3 . Written more explicitly, the kernel 
(151) is: 

^£2^3 = &(?2 +o- 3 -o-' 2 - a' 3 ) exp f ^- S(cr, cr')^ ; 



where 



wi(cr 2 - o"i) 2 + u 3 (cji - cj 2 ) 2 («iuiW2 + k 3 u 2 w 3 + KiK 3 niw 3 ) (fj 2 - a' x y 

S ( f7 ' a ) = =1 — } ; — \ 1 1 , — ^ 

u 2 x wi(k 1 ui +u 2 ) k 3 w 3 (k 1 u 1 + u 2 ) 

u 3 {k\u\w 2 + k 3 u 2 w 3 + K 1 K 3 uiw 3 )(a' 2 - a[) 2 , , 

H -, -7 r — 2(0-1 — a%)a 3 . (57) 

(K 3 U 3 W 3 + WlW 2 + U 3 W 2 )(KlUl + U 2 ) 



From Q53JI we get: 



, (wi + U 3 )w 2 (q2 - qi) + -U2-U3q3 + («2 + «lWl)u>lPl 

q 3 = 

^2^3 + ^2^1 + KiUltDi 
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and (|56|) becomes: 



dV <5(cx 2 + a 3 - 4 - 4) exp f^Efo a')) 4 £ V - O 



(wi + ii 3 )u 2 (4 - a' x ) + u 2 u 3 a' 3 + (u 2 + KiU^w^-^r 
<ro 5 (a -a) L x 

U 2 U 3 + U 2 Wi + KlUiWi 



x 5(4 + 4 - 4 - 4) exp ( ^£(4 a") ) . (58) 



Since 



(u 2 + kiUi)-^ = (w 1 + u s )u 2 a 1 -u 2 w 1 a 2 -u 2 u 3 a 2 + (wiu 2 + K 1 uiWi)a 3 , 

eq.(|58[) reduces to 

„ , „ ^3(^2 + ^3 ~ 4) + ("2^1 + Kimwi)4 „ „ 

o- 3 d(cr 2 + cr 3 - cj 2 - cr 3 J = ■ ■ d(a 2 + a 3 - a 2 - a 3 ). 

u 2 u 3 + U 2 Wi + KlUlWl 

□ 

Certainly, if this bosonic representation is used, the MTE involves integrations over R 
instead of summations over Zjv- The bosonic MTE may be proven directly with help of 
Gaussian integrations. 

There are indications that this bosonic model will not be critical all over and so be 
physically interesting [37] . 



6 Conclusions 

A large class of integrable 3-dimensional lattice spin models is constructed. The start- 
ing point and central object of the construction is an invertible canonical automorphism 
of a triple ultralocal Weyl algebra which satisfies two physically motivated principles: A 
current branching rule and a Baxter Z-invariance. To argue that this automorphism op- 
erator fulfils the tetrahedron equation requires no calculation. If the Weyl parameter q 
taken to be a root of unity, one chooses the N x N unitary representation of the Weyl 
operators and the N-th powers of the Weyl operators are represented by numbers. Accord- 
ingly, in this representation the automorphism operator can be decomposed into a matrix 
conjugation and a functional mapping. Quite trivially, the functional mapping operator 
satisfies a tetrahedron equation too. The conjugation matrix which can be expressed in 
compact fashion in terms of the Bazhanov-Baxter Fermat curve cyclic functions, satisfies 
the modified tetrahedron equation, (|28[). So we have a quantum problem with coefficients 
which are form a classical integrable system themselves. The large class of different so- 
lutions to the classical integrable system leads to a large class of quantum 3D-integrable 
models. The Zamolodchikov-Baxter-Bazhanov model is obtained if we choose the trivial 
classical solution. There is also a Gaussian continuum representation of the automorphism 
operator. 

Our expectation is that the framework set is broad enough to contain physically inter- 
esting 3D-integrable lattice models. An immediate application is to perform an asymmetric 
limit to construct new 2D-integrable chiral models involving rapidity parameters living on 
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managable algebraic curves. The vast possibilities to choose convenient parameterization 
of the 3D-Boltzmann weights should open new possibilities to solve the 2D-descendent 
models. 
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7 Appendix 

7.1 7^-mapping from eqs. (j^j) 

We give some details of the calculation leading from (|3* |) to We use the first four eqs. 
of © to express the currents |<^) [i = 1,2,3) and |<fo) in terms of \4>\) and |</> 3 ). This 
gives the four equations (not commuting any Weyl variables): 

\4> g ) ► (W 2 W1 -WjWa)^!) + (W1W3 - wiw 2 + g 1/2 W 2 U 3 - K3U 3 W 3 )|03) = 

|0e) > Q 1/2 (u 3 U 2 - U 2 U 3 )[(/)3) + |Kig _1/2 UiWi -U 2 Wi-g 1/2 U 3 (ui -U 2 )| |(/>i) = 

\<t>h) ► l( w 3 u l - u 'l w 3)(|^l) - |03>) 

+ | K 29 _1 ^ 2u 2 w 2 ~~ u 'l w 2 ~~ W 3 U 2 + l( w 3 u l + u 'l w 3)| {\4>l) + 103)) = 

\(f) f ) — > j/siuiwiwa - k 2 u 2 w 2 wi + g 1/2 K 3 u 3 w 3 (u 2 - Ul)| 

+ j^iuiw^w^ - w 3 ) + g 1/2 K3U 3 w 3 u 2 - g 1/2 K 2 U 2 W2U 3 | |0 3 ) = 0. (59) 

These homogenous eqs. allow the currents \(f>i) and |</> 3 ) to be chosen freely so that their 
coefficients must vanish and we have 

w 'l w 2 = w 2 w li U 3 U 2 = U 2 U 3! w 3 lu 'l = UiWj 1 (60) 

We see that for ultralocal Weyl algebras the three products wiw 2 , u 2 U3, U1W3 1 are 
centers of our mapping. From the second terms of \4> g ) and \<p e ) we get 

w[ = |w 2 wi - g 1 / 2 w 2 u 3 + k 3 u 3 w 3 ^ w 3 1 

u' 3 = (u 2 u 3 - g _1 / 2 u 2 wi + Kig _1 uiwi^ u^ 1 . (61) 

Using the center eas. (|60|) . we can express immediately also w 2 and u 2 in terms of the 
unprimed variables. 

In order to express in terms of the unprimed variables we take the difference of the two 
curly brackets in \cftf) and using ()60|) obtain 

/ / 1/2 / / / 1/2 



KiUjWjWa — q 1 k 3 u 3 u 1 w 3 = K2U2W2(wi — q 1 u 3 ) (62) 
Ising (|6T|) and u- = q 5i ' ] W- u- , u» Uj = q Si >i uj Uj leads to 

u[ |kiW 2 + QK 3 U 2 U ] ^ 1 W 3 + Q 1//2 KiK 3 W 3 | (Wi - g 1 / 2 ^) = K 2 U 2 W 2 (\Vi - ^ 2 U 3 ). (63) 
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One easily checks that (|6U|) . (|61|) . (|63|) can be written as (@J, ©• The inverse relations 
look similar, we give only one term: To express ui or W3 in terms of primed variables, we 
use the second term of \4>h), obtaining immediately 

ui = (uiw 2 + W3U2 - K2<?~ 1/2 u 2 w 2 ) w 3 _1 . (64) 
7.2 7^-mapping from "Linear Problem" 

It is often useful, instead of requiring (J3J), to start the derivation of the mapping 7^2,3 
postulating six "linear problem" -equations for "co-currents" (0j|. 

One may understand the nature of the mapping and the origin of the local Weyl algebra 
relations if we start from the most general consideration, i.e. introducing at each vertex 
three dynamic variables. 

In each sector 0j of the auxiliary plane there shall be a co-current (4>j\. At each vertex 
i some dynamic elements Uj, w$ and Vj acts on the four co-currents of the surrounding 
sectors such that 

1) the interaction with the co-current (0,-| of the sector between the arrows leads to a 
co-current (0j| (j ,1 ' 2 Uj at the vertex; with a co-current (0j| in the sector to the left of the 
arrows: to (0j| . With the (4>j\ of the sector to the right of the arrows the interaction 
produces (tfij | Vj , etc. 

2) The total co-current at a vertex is obtained by linear superposition of the four contri- 
butions and it must vanish, 

3) These properties are unchanged if one line is shifted through the vertex i (V — A- 
invariance) . 

The dynamic elements Uj, Wj and Vj form an associative noncommutative ring (actually a 
body). Co-currents belong to the formal left module of this ring. Note that the co-currents 
are not the conjugates of the currents used in eqs.(JHJ). 

Applied to the six vertices appearing in FigOl these rules lead to the following six equa- 
tions: 



(01 1 = 


(0c| 


+ ( 


:m <? 1/2 ui 


+ (<? 


>a| Wi 


+ 


(0e| Vi = 


0. 


(02 1 = 


(0a I 


+ ■ 


(0ek 1/2 U 2 


+ (<? 


^| w 2 


+ 


(0/| V 2 = 


0, 


(03 1 = 


(0c 1 


+ < 


!0ak 1/2 U3 


+ <? 


>&| w 3 


+ 


(03 1 V 3 = 


0, 


(01 1 = 


<0&| 


+ < 


:0,k 1/2 u; 


+ (<? 




+ 


(0/| Vi = 


0, 


(02 1 = 


(0c| 


+ < 


' / 1 1/2 / 

,0d| q 1 u 2 


+ (<? 


) 6 | w 2 


+ 


(0fc|v 2 = 


0. 


(03 1 = 




+ ■ 


(0ek 1/2 U 3 


+ (<? 


^1 w 3 


+ 


(0/1 v 3 = 


. 


is just a 


scale fact 


or and plays no 


role. 


The mapping 


is 



(65) 



demand of the linear equivalence of A-shape and V-shape triangles: any two equations of 
(|65[) must be linear combinations of the rest four. In the other words, any four external 
co-currents are to be linearly independent. 
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As an example of the application of this principle consider two linear combinations 
which are designed such that both do not contain the co-currents {4> a \, (4>c\, (4>h\ : 

(VI = (^-(^IKnY/V, 

(V = (Vslw^ 1 - (^w^ 1 + (0 2 |( Wl - g 1/2 u 3 )w 3 1 . (66) 
Explicitly, using (|65j) . we have 

(VI = {M + ^ a \<-^dWi\ l,2 ^i 

-(^IguX" 1 ^ + <<M(vi - v^wS-V^ui) , (67) 

= (4>b\ + (0 9 K V 3 W 3 1 + W 2 W 1 W 3 1 ~ ( 1 1/2 ™2 U 3 W 3 1 ) 

-(4>d\q 1/2 ^l^3 1 - (0e|(ViWg 1 + gu^W^ 1 - g 1 /2 U2WlW -l) 

+ (0/|v 2 (w 1 -g 1 / 2 u 3 )w3 1 . (68) 

The difference (VI — (V contains exactly four external co-currents, therefore due to the 
rank = 4 demand all the coefficients in (VI — must vanish identically. This gives in 
particular the expressions for and u 3 in the terms of the unprimed elements, and some 
other relations. 

In the same way one may consider all the other combinations of (<^| , . . . , (c/> 3 | giving 
linear combinations of four external co-currents. There are exactly eight independent 
relations for Uj , ... , v • . The solution of the linear equivalence problem may be written 
as follows: First, the six primed elements may be expressed explicitly 



where 



w 'l = w 2^3 J w 2 1 = w l * A 3 i U 3 V 3 1 = v l w l A 2 

w l v 'f 1 = V 3 U 3 lA 2 . u 2 _1 = u 3 lA l > u 3 = U 2 A 1 » 



Ai = u 3 u x 1 — q 1 ^ 2 w 1 u 1 1 + q 1 u 2 1 v 1 u 1 1 



(69) 



A.2 = u 3 v 3 1 w 2 v 2 1 + w-lV-l 1 u 2 v 2 1 — q 1 ^ 2 v 1 



A3 = w 1 w 3 1 — g 1 / 2 u 3 w 3 ~ 1 + w 2 ~ 1 v 3 w 3 _1 . (70) 
Also there exist three extra relations written below in the left column 

J „'-!„/ ..'-1 „-l„ „-l 



1 - u l w 3 vs U 3 V 3 v l w l = v l w l U 3 V 3 . 



vs 



/ /— 




U 3 V 3 




w^w 2 


= w 2 w 


U 3 U 2 


= U 2 U 3- 



1- 



V 3 W 3 1 V 2 W 2 1 = V 2 W 2 1 V 3 W 3 1 VS ii.Li,, lit li-,. 

(71) 

The right column in (|71[) is just an evident consequence of Q69JI . 

Two important notes are to be made. Recall, at the first, expressions (|^9*|) . (|70j) . (fTTl) are 
obtained in the framework of the associative noncommutative ring, i.e. we never changed 
the order of the invertible elements Uj, Wj, Vj, u£, w- and v^. Secondly, since the linear 
equivalence principle implies eight independent relations, there are eight independent re- 
lations between nine relations ()69j) . (|71[) . In the framework of the noncommutative ring 
the proof of this miracle is a quite complicated exercise. 
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Now the question is how to get rid of one degree of freedom in the primed elements. 
The answer lies in the comparison of the left and right columns of (|71|). One may see that 
it is natural to identify left and right columns. The identification produces two demands 

v^ur 1 ~ w< , v^r 1 ~ Ui (72) 

and the analogous relations for the primed elements. Evidently, eqs. (|72[) fix the Weyl 
algebra structure. 

Let now the initial elements Uj, Wj and v, obey the local Weyl algebra relations © with 
Vj = KjUjWj, Ki G C. Then (|69j ) .l|71 |) immediately give v£ = k^u^w^, where k\ are centers 
such that and the mapping (|69j) defines the automorphism of 

the local Weyl algebra. External "physical" demand is that since Ki are centers, a well 
defined automorphism should not change them: k[ = K i . After the identification of Ki and 
k\, eqs. (|p^|) . l(7n|) coincide with ©,© up to the rescaling A2 = q , kiK3A2. 
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